Structure formation in a nonlocally modified gravity model 
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We study a nonlocally modified gravity model proposed by Deser and Woodard 1 which gives an 
explanation for current cosmic acceleration. By deriving and solving the equations governing the 
evolution of the structure in the Universe, we show that this model predicts a pattern of growth that 
differs from standard general relativity (+dark energy) at the 10-30% level. These differences will 
be easily probed by the next generation of galaxy surveys, so the model should be tested shortly. 



PACS numbers: 

I. INTRODUCTION 

There are multiple threads of evidence - from distant 
supernovae [2J [3] , the cosmic microwave background [U- 
7 , the galaxy distribution [8], and correlations of these 
with each other [5] - that the Universe is accelerating 
The multiple sets of evidence argue against any one sys- 
tematic affecting the qualitative conclusion that a > 
where a is the scale factor of the Universe. 

General relativity maintains intact our Newtonian in- 
tuition that gravity should lead to deceleration, in the 
case of the smooth Universe that a should be negative. 
The one caveat is that general relativity allows not only 
energy, but also pressure, to impact geometry, and a sub- 
stance with negative pressure can drive acceleration. One 
possibility then is that the multiple strands of cosmo- 
logical evidence point to the need for a new, previously 
undiscovered substance with negative pressure, so-called 
dark energy [10H15] . 

Another possibility is that general relativity is 
wrong [16) . The theory was developed using informa- 
tion and intuition honed on relatively small scales (the 
Solar System), so it would not be shocking if it needs 
to be modified on the larger scales probed by cosmology. 
There have several problems so far with this approach. 
Perhaps, most generally, we do not have a set of simple 
guidelines, thought experiments, and arguments of ele- 
gance that motivated Einstein as he constructed general 
relativity. Even if beauty is abandoned, a modification of 
gravity must still confront three problems, two of which 
are related. First, almost all models contain a mass scale 
that must be set to be much smaller than any mass found 
in nature, less than 10 -33 eV. Why such a small mass 
scale exists, and how it can be protected in the presence 
of even weak interactions with the rest of physics, is a 
profound problem. Related to this is the fine tuning prob- 
lem in time that afflicts these models: the modifications 
to gravity just happen to be important today, not at any 
time in the past. Believing that we live in a special time is 
yet another price we pay for these modifications. A third 
problem with modified gravity models is that they must 
be introduced in a way that does not violate the successes 



of general relativity in the Solar System. Indeed, these 
constraints |17j doomed one of the first modified gravity 
models introduced to explain acceleration and still place 
tight constraints on many others (e.g., see [T8H2U] ). 

One class of modified gravity models that rises above 
these challenges contains non-local interactions p], i.e., 
terms in the Lagrangian that depend on the values of 
the fields at more than one point in space-time. In par- 
ticular, Deser and Woodard PQ pointed out the utility 
of terms that are functionals of \3~ 1 R, where □ is the 
d'Alembertian and R the Ricci scalar. In the cosmolog- 
ical context, C~ 1 R grows very slowly: as {t/teq) 1 ^ 2 in 
the radiation dominated era and logarithmically in the 
matter dominated era. So, at the time of Nucleosynthe- 
sis n _1 i? is about 10 -6 and at matter-radiation equal- 
ity it is only order one. These terms therefore naturally 
are irrelevant at early times and begin to affect the dy- 
namics of the Universe only after the matter-radiation 
transition, thereby circumventing some of the fine tuning 
difficulties. Since is dimensionless, the functional 

that multiplies R has no new mass parameter. Finally, 
because n _1 i? is extremely small in the Solar System, 
these models easily pass local tests of gravity. 

There is some theoretical motivation for these kinds of 
terms. Away from the critical dimension in string theory, 
R\3~ 1 R is precisely the term generated by quantum cor- 
rections, as first pointed out by Polyakov 21 . More gen- 
erally, there have been several arguments for considering 
non-local theories [32H2S] • It might be possible to rewrite 
these models in terms of local models with an auxiliary 
scalar field [57) [58]. Given that they can explain the 
current epoch of acceleration, it is not surprising that 
an early epoch of inflation has also been attributed to 
non-local interactions [29H3T) . A realistic model for ac- 
celeration was constructed in Ref. [32J, where they chose 
the arbitrary function of so that it reproduced the 

exact expansion history generated by standard ACDM. 
This working model demonstrated the feasibility of the 
non-local approach to generating acceleration and there- 
fore we focus on it here, but of course it also points to 
a problem inherent in this way (and many others) of 
modifying gravity: there is a completely arbitrary free 
function that can be chosen at will to fit the expansion 
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history. 

One of the most intriguing aspects of modified gravity 
models in general is that, even if they arc constructed to 
reproduce the expansion history equivalent to that of a 
given dark energy model (such as ACDM) , perturbations 
will often evolve differently than in a model with stan- 
dard general relativity and dark energy. This has led to 
the idea that the way to distinguish dark energy models 
from modified gravity models is to measure the growth 
of structure in the Universe. While there has been some 
work to date on perturbations in non-local models |33j . 
there have been no concrete predictions for the growth of 
structure and how this may differ from dark energy mod- 
els. Here we work out the perturbation equations in this 
new class of models and solve perturbatively to arrive at 
a concrete prediction for two observables as a function 
of redshift. The deviations from dark energy models are 
at the 10-30% level and have a characteristic signature 
as a function of redshift, which suggests that the class of 
models will be tested by upcoming surveys. 

The rest of the paper is organized as follows: SjTT] re- 
views the model of Ref. [32] and presents the equations 
governing the expansion history in this model. The next 
section perturbs the model around this zero order solu- 
tion to arrive at the equations that govern the growth 
of structure. The next two sections solve these equa- 
tions perturbatively for two quantities: $ + $ and G c ff 
where and $ are the two scalar potentials in the per- 
turbed Friedman Robertson Walker metric and G c ff is 
the proportionality constant between V 2 $ and 4irp5 in 
the Poisson equation. In general relativity, >F + $ = 
and G c ff = Gn, Newton's constant, but in modified grav- 
ity models these equalities often do not hold. Indeed, we 
show in SjV] and §IV| that the deviations in this non-local 
model are significant and measurable. 



II. MODEL AND ZEROTH ORDER EQUATION 

We take a model in which the Einstein-Hilbert La- 
grangian is nonlocally modified as pQ 



16ttG 



(1) 



where / is a function of its dimensionless argument that 
will be determined by matching the expansion history of 
ACDM. To understand how this new term helps drive 
acceleration today, note that in a smooth expanding 
Universe described by the Friedmann-Robertson- Walker 
metric, the non-local operator at time t is a double inte- 
gral over all previous times: 

Xsn " B =-/'^oi'' i<v < t "' s ( t "' 

* fit' t*' r . 

I dt"a 3 (t")\6H(t") + 12H 2 (t'"> 



a 3 (f) 



where here and throughout we denote quantities in the 
smooth FRW background with over-bars, e.g. R, to dis- 
tinguish them from the full perturbed quantities; and 
H = a/a is the expansion rate. In a perfect radiation 
dominated Universe, R — 0. So, in the standard cos- 
mological picture, with radiation giving way to matter 

at a EQ = 3 x 10~ 4 , X = □ ~ R is very small until agQ 
after which it grows logarithmically so remains of order 
unity today (solid curve in Fig. [I]) . During the epoch in 
which the Universe is purely matter dominated, a cx t 2 ^ 3 

so n -R| MD = — 21n(a/aEQ), an approximation that 
Fig. [T] shows captures the growth at late times. 
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FIG. 1: The evolution of X = C R. Dashed line shows the 
analytic solution in the matter era: X — — 21n(a/aEQ). The 
flat (green) curve shows the function f(X) chosen in Ref. |32| 
to fit the ACDM expansion history. Its amplitude and shape 
are chosen so that it has negligible impact on the expansion 
until relatively recently. 

It is straightforward to choose a function of X = n _1 i? 
in the action of Eq. (JlJ that begins to modify the expan- 
sion history only at late times when its argument becomes 
of order unity. To see what form the function must take, 
we vary the action with respect to the metric, obtaining 
the modified Einstein equations: 

G„„ + AG M „ = 8nGT^ . (3) 

Here the correction to the Einstein tensor, AG M „ is pQ 



AG™ 



Gfj, u + g^ v LJ 




Applying these general equations to the FRW metric 
leads to 



(2) 



-2H - 3H 2 + 



3H 2 
1 
3^2 



+ AG 00 
5 ij AGn 



8irGp 
SttGP. 



(5) 
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Here p and P are the energy density and pressure and 
the zero order (i.e, homogeneous) nonlocal terms are pQ 



AG 00 = 



3H 2 + 3Hd t ] | / (x) += [nf (x 
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Carrying out the reconstruction program using this re- 
vised Friedmann equation, Deffayet and Woodard |3"2"] 
showed that the function 

f(X) = 0.245 [tanh(0.350y + 0.032F 2 + 0.003F 3 ) - 1 

(9) 

where Y = X + 16.5 leads to an expansion history equiv- 
alent to that produced by ACDM when the rest of the 
cosmological parameters are set to their default values, 
{f2curv,fim,fir} = {0,0.28,8.5 x 10~ 5 }. This function, 
called the non-local distortion function, is depicted in 
Fig. [I] which gives a sense of the fine-tuning needed in 
this model. The function must be very flat so that the 
new terms do not alter the expansion history until only 
recently. The time at which the modifications kick in and 
the amplitude are both fine-tuned to produce the correct 
expansion history. Due to the logarithmic nature of the 

growth of n R, though, this fine-tuning is many orders 
of magnitude less severe than that in most other models. 
Apart from aesthetic arguments, the model defined by 
this function is concrete enough to enable predictions to 
be made for the growth of structure. 



In our notation, little 5 refers to first order perturba- 
tions, and capital A means the correction to the original 
Einstein equation due to the non-local terms. The time- 
time component leads to the modified Poisson equation 
on small scales, and we focus on that modification first. 
The space-space components can be arranged to yield an 
equation relating the two scalar potentials, and we treat 
the non-local modification of that equation afterwards. 



A. Modified Poisson Equation 

Some of the terms in the time-time component of 
Eq. (11) arc standard; for example, in Fourier space, 



SG oo = 
ST o = 



6H$ 
pS. 
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On the first line, we drop the first term because we are 
working in the sub-horizon limit so the second term, 
which is of order (k/H) 2 with respect to the first, domi- 
nates. On the second line, p is the mean density of mat- 
ter, and S the over-density, (p — p) /p. With just these 
terms, we recover the standard Poisson equation. 

The new terms arising from the non-local distortion 
function, SAG, however modify the Poisson equation. To 
determine these new terms, we see from Eq. Q that we 
need to evaluate the perturbed AGqo or 
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where 



III. PERTURBATION EQUATIONS 

In order to follow the growth of inhomogeneities in the 
universe, we write the metric as the background FRW 
geometry plus small perturbations. Plugging the metric 
into the field equation (|3| and expanding it to the first 
order gives the evolution equations for the perturbations. 
We consider only scalar perturbations and follow the no- 
tation of Ref. pS] (Eq. (4.9) there): 

goo(t,x) = -l-2$!(t,x) , g Oi (t,x) = 0, 

g l3 {t,x) = ^-a 2 (i)[l + 2$(i,f)] . (10) 

With these ingredients the modified field equations at 
first order are 



AG 00 
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(11) 
(12) 



(15) 



Capturing all of the relevant first order terms is, for 
the most part, standard but quite tedious. The inverse 
d'Alembertian's however introduce a bit of a twist, so we 
will evaluate the first term as an illustration and then 
present the final equation. 
The first term is 



SAGn 



5 (Poo/) = V Q0 6f + SV QO f(X) . (16) 



first term 



The simplest piece of this is the one depicted in Fig. [T] 
f(X), which - using Eq. ^ and Eq. Q9j) - is a given 
function of time. The zero order derivative operator is 
also straightforward; plugging in the FRW metric leads 
to 



V 



oo 



3H 2 + ZHd t 
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(17) 
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where the first two terms have been dropped because 
they are smaller than the third by a factor of (H/k) 2 
(the sub-horizon limit again). Slightly more complicated 
is the first order perturbed derivative operator 



k 2 V 2 
SV 00 = 6i7$ + 3$a i + 2- ¥ $ + 2($-*)- 

a 1 



i$ 



kkd k 



(18) 



where the sub-horizon limit kills the first two terms, and 
the two terms containing spatial derivatives have been 
set to zero because they vanish when operating on the 
homogeneous f(X). We therefore have 



SAG, 



mi 



= %f'{X)SX + 2%<S>f(X) (19) 



first term 



where we've set Sf = f X SX. 

To obtain SX, we expand both the d'Alembertian and 
the Ricci scalar, so that 



n^i? = n r + h sr-h (sn)x. (20) 

The last two terms on the right constitute SX, so that 



sx = n [sr - (sn)x\ 



(21) 



Both SR and SO can be written in terms of the scalar 
perturbations. We find in the sub- horizon limit 



SR 

sn = 0. 



v 2 

-2-=-(tf + 2$) 



(22) 



To evaluate □ SR, it is useful to introduce the retarded 
Green's function that satisfies 

UG iet {x;x') = ( K -d 2 -Wd t + ^G Iet (x;x') = 5\x-x') 

(23) 

Once we solve for the Green's function, SX in the sub- 
horizon limit can be written as 



SX(x)= / <Fx' G ret {x,x')5R{x'). 



(24) 



In general the Green's function can be constructed us- 
ing the massless, minimally coupled scalar mode func- 
tions u(t,k) for arbitrary a(t): 

„ , 0{t-t')a 3 {t') f d 3 k $.rg_2) 

u(t, k)u* (t' } k) - u* (t, k)u{t', k) (25) 

Although there is no general solution for G rot , the Green's 
function can be captured analytically in the small scale 



limit (in Fourier space, k/Ho —} oo) . Fortunately, this 
is precisely the limit we are interested in, since it is these 
scales that can be most easily probed by observations. In 
this limit, using the WKB approximation we find 



u(t,k) 



1 exp 



2k 



a(t) 



(26) 



so that the integral over d x in Eq. (24) leads to an 



expression in Fourier space for the perturbation to X: 



SX(k,t)= / dt'G{t,t';k)SR(k,t') 
Jo 

where the Fourier space Green's function is 



G{t,t';k) 



1 



ka(t) 



a 2 (t') sin 



dt" 



v a(t") 



(27) 



(28) 



Using Eq. (|22|) leads to 
SX(k,t) 



k 2 f* 

2— / dt'G{t,t';k) ^{k,t') + 2^{k,t') 

(29) 

Treating the other terms in Eq. (|14j) similarly leads to 



our final expression for the modified Poisson equation in 
the sub- horizon limit: 

fc 2 $ + k 2 E[<$>] = 4:TTGa 2 pS. (30) 

where the new terms are functionals of and 



£?[*,*] = $J/(X) + g 



Rf'iX 



l -\l'{X)sx + n 1 j'[x)sn } c-jd 



Note that all the new terms are of order fc 2 $ since the 
inverse d'Alembertian scales as k~ 2 when acting on the 
spatially varying SR but as H~ 2 when acting on the ho- 
mogeneous Rf'(X). 



B. Gravitational Slip Equation 

The equation for the so-called gravitational slip (351 — 
157] . the difference between the amplitudes of the two 
scalar potentials follows from the space-space compo- 
nents of Einstein's (now modified) equations. To derive 



1 In the limit that C acts on a spatially homogeneous field 
(in Fourier space, k/H — > 0), the Green's function is 5 3 (x — 

x') f* dt" a JrnT ■ Plugging this expression into the equivalent 

of Eq. for _R and changing the order of integration leads 

directly to Eq. |2l. 



5 



the modified equation for the gravitational slip, we start 
with the p ertu rbed space-space components of Einstein's 
equation, (12 1. Contracting with kk? — (l/3)(5 lJ to ex- 



tract the longitudinal traceless component leads to 



3a 2 



k 2 ($ + *)- 



SAG t] = 8ttG (p - 



k l W - (1/3)<T J 

(32) 

where a (in the notation of Ref. [35]) represents the 
anisotropic stress. At late times, when contributions to 
the energy budget from relativistic species such as pho- 
tons and neutrinos are negligible, the anisotropic stress 
vanishes, so $ + ^ = is a robust prediction from general 
relativity. 

To quantify the extent to which this prediction is vio- 
lated in the nonlocal model, we must compute the new 



P)C7 



terms in Eq. (32), those proportional to SAGij. For the 



most part, evaluating these requires the same set of tech- 
niques used in arriving at the modified Poisson equation. 
Some of the new ingredients (as always in the sub-horizon 
limit) are: 



— j k ~\~ k<i k j 
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*) +kikj (* + *). (33) 



Dropping the anisotropic stress, we arrive at the modified 
equation for the gravitational slip: 



($ + *) 



-(* + ¥)J/(jr) + i= 
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-f(X)SX- =f(XjSR 



(34) 



IV. SOLUTION OF THE POISSON EQUATION 

To solve the modified Poisson equation, note that the 
new terms due to the non-local interactions are small. 
This follows from the amplitude of / as depicted in Fig. [j] 
(note that the plotted curve is 10 x /). We can therefore 
solve the equation perturbatively by inserting the stan- 
dard solution for $ and into the new terms and deter- 
mining their impact on the Poisson equation. This ap- 
proximation works as long as the effect of the new terms 
remains small. 



We can write Eq. (30) as 



fc 2 $ = 4irG e ga 2 p5 
so that the fractional change G e g/G is 



G c ff 
~G 
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(36) 



The changes to Newton's constant are therefore cap- 
tured by evaluating the integrals that define E[$] in 
Eq. pill. Using Eq. d22|, Eq. p9| and setting * = 



in accordance with our strategy of perturbing around the 
GR solution leads to 



xG{t,t';k) 



Rf(X 
$(M' 



k- 



dt' 



o a 2 (*0 
[f'(X(t)) + f'(X(t'))] 



(37) 



To evaluate the integral we use standard linear time de- 
pendence of the gravitational potential, as encoded in the 
growth function: 



$(£;,a) = <f>(k,a = 1) 



Dx(a) 



aDi (a = 

where the growth function in ACDM is 

da 1 



Di(a) 



a' 3 i?(a') 3 ' 



(38) 



(39) 



In Figure [2] two cases k = lOOiJo and k = lOOOiJo are 
plotted. There is little fc-dependence, at least on these 
large scales that are close to linear. The deviations from 
GR though are quite apparent at late times and should 
be detectable by upcoming surveys. 




FIG. 2: The fractional change to Newton's constant, —gp- as a 
function of redshift z. The two curves, which depict the evo- 
lution for k — 0.03, 0.3h Mpc -1 , show that the modification 
is virtually scale-independent, at least in the linear regime. 



V. SOLUTION FOR THE GRAVITATIONAL 
SLIP 

In this section we calculate the gravitational slip ($ + 
^) in the nonlocal model, starting from Eq. (34) and 



inserting the GR solution for the potentials in all the 
new terms on the right hand side. The new terms then 
become 



* + $ 
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dt'G(t,t';k) 
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FIG. 3: Gravitational slip as a function of redshift in the 
nonlocal model. The two curves, barely distinguishable, are 
for k = 0.03 h Mpc" 1 and it = 0.3 h Mpc" 1 . 

Fig. [3] shows the evolution of the gravitational slip as 
a function of redshift. Here too we find a detectable 



deviation from the prediction of GR. 



VI. DISCUSSION 

A nonlocal model of gravity is potentially interesting 
in that it can explain the observed acceleration of the 
Universe without excessive fine-tuning. In this paper, 
we have showed that it makes predictions for the growth 
of structure that differ from those of models based on 
general relativity. These differences are at the level that 
the model should be tested by upcoming galaxy surveys. 



Acknowledgments 

We are grateful for very helpful discussions with 
Matthew Dodelson, Mark Trodden, and Richard 
Woodard. The computation of the integrals was car- 
ried out with the aid of Mathematica. This work was 
supported by the U.S. Department of Energy, including 
grant DEFG02-95ER40896, and by the National Science 
Foundation under Grant AST-090872. SP was supported 
by a Fermilab Fellowship in Theoretical Physics. 



[16] 



S. Deser and R. Woodard, Phys.Rev.Lett. 99, 111301 [17 
(2007), 0706.2151. [18 

A. G. Riess et al., Astron.J. 116, 1009 (1998), astro- 
ph/9805201. [19 
S. Perlmutter et al., Astrophys.J. 517, 565 (1999), astro- 
ph/9812133. [20 
S. Dunkley et al., Astrophys.J. Suppl. 180, 306 (2009), 
0803.0586. [21 
E. Komatsu et al., Astrophys.J. Suppl. 180, 330 (2009), [22 
0803.0547. [23 

B. D. Sherwin, J. Dunkley, S. Das, J. W. Appel, [24 
J. R. Bond, et al, Phys.Rev.Lett. 107, 021302 (2011), 
1105.0419. [25 
A. van Engelen, R. Keisler, O. Zahn, K. Aird, B. Benson, 

et al., Astrophys.J. 756, 142 (2012), 1202.0546. [26 

A. G. Sanchez, C. Scoccola, A. Ross, W. Percival, [27 
M. Manera, et al. (2012), 1203.6616. 

R. Scranton et al. (SDSS Collaboration) (2003), astro- [28 
ph/0307335. [29 
J. A. Frieman, C. T. Hill, A. Stebbins, and I. Waga, 
Phys.Rev.Lett. 75, 2077 (1995), astro-ph/9505060. [30 
K. Coble, S. Dodelson, and J. A. Frieman, Phys.Rev. 
D55, 1851 (1997), astro-ph/9608122. [31 
R. Caldwell, R. Dave, and P. J. Steinhardt, [32 
Phys.Rev.Lett. 80, 1582 (1998), astro-ph/9708069. 

B. Ratra and P. J. E. Peebles, Phys.Rev. D37, 3406 [33 

(1998) . 

C. Wetterich, Nucl. Phys. B302, 668 (1998). [34 

D. Huterer and M. S. Turner, Phys.Rev. D60, 081301 

(1999) , astro-ph/9808133. [35 
S. M. Carroll, V. Duvvuri, M. Trodden, and M. S. Turner, 
Phys.Rev. D70, 043528 (2004), astro-ph/0306438. [36 



T. Chiba, Phys.Lett. B575, 1 (2003), astro-ph/0307338. 
S. Nojiri and O. S. D., Int. J. Goem. Meth. Mod. Phys. 

4, 115 (2007). 

T. D. Saini, S. Raychaudhury, V. Sahni, and A. Starobin- 
ski, Phys. Rev. Lett. 85, 1162 (2000). 
Z. K. Guo, N. Ohta, and Y. Zhang, Phys. Rev. D72, 
023504 (2005). 

A. M. Polyakov, Phys.Lett. B103, 207 (1981). 

L. Parker and D. Toms, Phys.Rev. D32, 1409 (1985). 

T. Banks, Nucl.Phys. B309, 493 (1988). 

H. W. Hamber and R. M. Williams, Phys.Rev. D72, 

044026 (2005), hep-th/0507017. 

T. Biswas, A. Mazumdar, and W. Siegel, JCAP 0603, 
009 (2006). 

N. Barnaby, Nucl.Phys. B845, 1 (2011), 1005.2945. 

5. Capozziello, E. Elizalde, S. Nojiri, and S. D. Odintsov, 
Phys.Lett. B671, 193 (2009), 0809.1535. 

N. Koshelev, Grav.Cosmol. 15, 220 (2009), 0809.4927. 
N. Tsamis and R. Woodard, Annals Phys. 267, 145 
(1998), hep-ph/9712331. 

S. Nojiri and S. D. Odintsov, Phys.Lett. B659, 821 

(2008), 0708.0924. 

T. Prokopec (2011), 1105.0078. 

C. Deffayet and R. Woodard, JCAP 0908, 023 (2009), 
0904.0961. 

T. S. Koivisto, Phys.Rev. D78, 123505 (2008), 
0807.3778. 

S. Dodelson, Modern Cosmology (Academic Press, San 
Diego, 2003). 

S. F. Daniel, R. R. Caldwell, A. Cooray, and A. Mel- 

chiorri, Phys.Rev. D77, 103513 (2008), 0802.1068. 

S. F. Daniel, R. R. Caldwell, A. Cooray, P. Serra, A. Mel- 



7 



chiorri, et al., Phys.Rev. D80, 023532 (2009), 0901.0919. [38] C.-P. Ma and E. Bertschinger, Astrophys.J. 455, 7 
[37] E. Bertschinger, Phil.Trans.Roy.Soc.Lond. A369, 4947 (1995), astro-ph/9506072. 

(2011), 1111.4659. 



